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ABSTRACT 



This paper deals with black holes, bubbles and orbifolds in Gauss- 
Bonnet theory in five dimensional anti de Sitter space. In particular, we 
study stable, unstable and metastable phases of black holes from ther- 
modynamical perspective. By comparing bubble and orbifold geometries, 
we analyse associated instabilities. Assuming AdS/CFT correspondence, 
we discuss the effects of this higher derivative bulk coupling on a specific 
matrix model near the critical points of the boundary gauge theory at 
finite temperature. Finally, we propose another phenomenological model 
on the boundary which mimics various phases of the bulk space-time. 
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1 Introduction 



Black hole in AdS space has a remarkable property that it undergoes Hawking-Page 
(HP) phase transition [1]. Asymptotically AdS space allows two kinds of black hole 
configurations. By comparing their sizes with respect to the AdS scale, one can 
characterise these holes. While the small black holes have horizon sizes less than 
that of the AdS scale, the big black holes are larger than the AdS scale. These small 
black holes, however, are unstable with negative specific heat; leaving big black holes 
as the stable configurations in AdS space. Furthermore, it was noticed that, as we 
tune the temperature close to the inverse of the AdS scale, there is a first order phase 
transition. At a temperature below the inverse AdS scale, system prefers thermal AdS 
space, while at higher temperature, it is the big black hole phase which minimises 
the energy of the system. This crossover is known as HP transition. Via AdS/CFT 
correspondence [2], this phenomenon was found to have its imprint on the gauge 
theory residing on the boundary. Witten argued [3] that, on the boundary, the HP 
transition represents a large N deconfinement transition of the gauge theory at strong 
coupling. 

Certain analytical continuation of the black hole metric in AdS space gives bubble 
of nothing solution [4,5]. These are the analogues of Witten's Kaluza-Klein bubbles 
in flat space-time [6]. Bubble spacetime corresponds to time dependent configura- 
tion and, as we review later, in five dimensions, the boundary metric is dS 3 x S 1 . 
By using AdS/CFT correspondence, one then hopes to learn about gauge theory on 
time-dependent geometries. An important ingredient in the AdS/CFT correspon- 
dence is the principle of holography [7] . According to this principle, the physics of a 
gravitational theory is dual to a different theory in one lower dimension. Conversely, 
given a dual theory on a boundary, we must consider all possible bulk space-times 
whose boundaries have the same intrinsic geometry as the background of the dual 
theory. Now, given that the boundary metric is dS^ x S 1 , there exists another bulk 
geometry known as AdS orbifold. These are the five dimensional analogues of BTZ 
black holes [8,9]. These orbifolds are however unstable and below a critical size of the 
boundary S 1 , they decay to bubbles of nothing [10,11]. A bubble, once formed, ex- 
pands exponentially and fills up the whole space-time. Though in AdS, this is similar 
to the decay of the Kaluza-Klein space to nothing. For other studies involving de- 
cay into bubbles of nothing, obtained by analytically continuing black hole solutions 
see [12-14]. 

In this paper, after briefly reviewing the HP transition and orbifold decay in 
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section 2, we analyse the response of these phenomena as we perturbatively increase 
the gravitational strength . Our study is partly motivated by recent works in [15-17]. 
In these papers, authors have argued in different ways that a version of HP transition 
occurs even at weak coupling gauge theory. By AdS/CFT dictionary, this would 
show up as a transition in strongly coupled gravity theory in the bulk. Noting the 
fact that string theory in AdS space is as yet poorly understood, we study a much 
simpler system in this paper. We add higher derivative terms in the supergravity 
action and study their effects on HP transition as well as on orbifold decay. We note 
here that higher derivative terms would arise in gravity action due to a 1 corrections 
in underlying string theory. While a study with a general class of higher derivative 
terms would be desirable, in this work, we consider only the effects due to Gauss- 
Bonnet (GB) terms. One advantage of working with GB correction to the gravity 
action is that the black holes, bubbles and orbifolds can be constructed explicitly. 

In section 3, we analyse the black holes in GB theory with a particular focus on 
their phase structures in five space-time dimensions. The phase structure depends 
crucially on the GB coupling. For certain range of coupling, there exists three black 
hole phases. We call them, small, intermediate or unstable and big black hole phase. 
It turns out that there are two first order phase transitions. One of them is from small 
black hole to the big one at a temperature scale much lower than that of inverse AdS 
curvature. The other one is similar to that of usual HP transition where a crossover 
occurs from thermal AdS to the big black hole phase. We compute the change in HP 
temperature in powers of the GB coupling at the crossover. 

In section 4, we study the bubbles in GB theory. We find that there exist bubble 
of nothing solutions for any value of the asymptotic circle. This is unlike the case 
in simple AdS gravity, where the bubble to exist, the circle size needs to be less 
than a critical value. Consequently, by computing energy densities of the bubble and 
orbifold spacetimes, we argue that orbifolds are always unstable and decay to bubble 
of nothing. The decay rate can then be easily computed by identifying the bounce 
solution. 

Many papers in the recent past have analysed the partition function of free M = 4 
super Yang-Mills theory and argued that the large N deconfinement transition occurs 
even at zero coupling [16, 17]. In fact, it turned out that the transition appears 
exactly at the Hagedorn temperature of the low temperature thermal AdS phase. 
Subsequently, non-perturbative 1/N effects near the Hagedorn transition was studied 
in [18]. This has been analysed in the other limit of the 't Hooft coupling, A — > oo by 
proposing a phenomenological matrix model [19,42,43]. 
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It is known that GB corrections arise in heterotic string theory, see for example 
[20]. In type IIB theory the corrections start off with R A . In the case with Gauss 
Bonnet terms we do not expect the boundary theory on S 3 x S* 1 to be Af = 4 Yang 
Mills theory. However in the limit a — > the boundary theory should reduce to 
the strongly coupled SYM theory. With the R 2 corrections turned on, the gravity 
theory should correspond to some deformation of Af = 4 SYM. In Section 5, we study 
this effective theory by using a phenomenological matrix model proposed in [19]. 
This model is characterised by two parameters which we call, following [19], a and b. 
Generally, (a, b) depend on the gauge theory temperature and the 't Hooft coupling 
A. Following AdS/CFT, the effect of adding higher derivative terms in the bulk 
translates to A corrections to the boundary gauge theory. Assuming an universal 
nature of the (a, b) model around the critical points, we analyse the A dependence of 
parameters (a, b) around the HP points. We do this numerically in section 5. 

Finally, in section 6, we construct a toy model which captures the whole phase 
diagram of the bulk. However, this requires introduction of four parameters in the 
matrix model potential. These four parameters again depend on the temperature as 
well as the gauge coupling. We then study the qualitative behaviour of this model. 
This paper ends with a discussion of our results. We hope to report on a similar 
analysis for the type IIB theory with i? 4 term in a future publication [21]. 

2 Black hole, bubble and AdS orbifold 

In this section, we briefly recall black holes, bubbles and AdS orbifolds in AdS grav- 
ity in five dimensions. We review the instability associated with the AdS orbifolds 
and calculate the decay rate of these orbifolds to bubbles in the supergravity limit. 
Furthermore, we compute an enhancement of the decay rate due to string wrapping 
around a compact direction of the AdS orbifold. We also briefly review, following a 
suggestion due to Horowitz [22], as to how this decay may be catalysed by tachyon 
condensation in string theory. 

Black Hole: 

Black hole in pure AdS-gravity is parametrised by a single parameter associated 
with the energy of the hole. Denoting this parameter as to, we may write the metric 
of the black hole as 

ds 2 = -(1 + £ - ™) dt 2 + (1 + £ - + r 2 (d6 2 + coste 2 ,), (1) 

L T' I T' 
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where / is related to the cosmological constant present in the action. This metric 
asymptotically approaches AdS space. The singularity is at r = and the horizon is 
located at r+, where r + is a solution of equation 

The Euclidean version of the metric is free of any conical singularity if the Euclidean 
time has certain periodicity. Defining 

near r = r + , we can write the metric as 

9 /r+ m\2 
i -t- i — — 



ds = d P + (iT + ^ j P d X + rldnl (4) 

where, we have used t = %x- From here it follows that the metric is conically non- 
singular if x has a period 

2nl 2 r 3 + 2nr + l 2 

Inverse of this periodicity is then identified with the black hole temperature. We note 
that, for a fixed temperature (above a certain critical value), we always get two black 
hole solutions with two different horizon sizes. We distinguish these two by calling 
them small and big black holes. At the critical temperature, both these black holes 
meet. It turns out that the smaller black hole is unstable with negative specific heat 
while the bigger one is stable. The free energy of the black holes is given by, 

0-2 2 2 

Here, k 5 is related to the five dimensional gravitational constant. We note that if the 
size of the black hole is larger than the AdS scale /, the free energy becomes negative. 
Since this is less than the free energy of thermal AdS space, there is a first order phase 
transition. From (5), we see that for r + = I, T = T c = So the transition occurs 
from thermal AdS space phase to the black hole phase as we increase the temperature 
beyond T c . The crossover between these two geometries is known as HP transition [1]. 
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Bubble: 



The other space-time that is of our interest is the AdS bubble [4,5, 10, 11]. The 
metric can be obtained by analytically continuing ( t — > %x and 6 — > n/2 + it) the 
black hole solution given in (1). We get 

ds 2 = (1 + £ - ™ 2 )d X 2 + (1 + £ - ^)-Mr 2 + r 2 (-rfr 2 + coshVdQ 2 ,). (7) 

If x is restricted to the period as in (5), the metric is non-singular for r > r + . This 
geometry is known as a bubble of nothing solution. For large r, at any time r, the 
metric is x circle times a two sphere. This circle collapses at r = r+. However, the 
two sphere approaches a finite size r+coshV. This two sphere is the boundary of the 
bubble. We see that the metric is time dependent and is asymptotically dS 3 x S 1 . 
As can be seen from the dotted line in Figure 6, below a certain critical value of Ax, 
for a given Ax, there are two possible bubble solutions. The smaller one, however, 
is expected to be unstable as its Euclidean continuation suffers from having modes 
with negative mass 2 [24]. The critical size of the x circle, above which there are no 
bubble solutions, is given by 

AXc = 7^ forrc = 7f (8) 

AdS orbifold: 

The AdS orbifold, that we consider here, has been discussed in great detail in 
[8,9,23]. A five dimensional AdS space is defined as the universal covering space of a 
surface which obeys 

2i 2i 2i 2i 2 2 i2 (q\ 

where, as before, / is the AdS curvature radius. The orbifold is obtained by simply 
identifying points along the boost 

£ = r -±( Xi d 5 + x 5 d 4 ), (10) 

where r + being an arbitrary constant. Since the norm of the boost is given by 
£ 2 = r \{— x \ + x \)/l 2 , £ 2 can be positive or negative. However, to avoid closed 
time-like curves, the region £ 2 < is removed from space-time. In an appropriate 
coordinate system this orbifolded space can be represented as 

rfs 2 = ( ~ 2 _ r 2 + ilcosh 2 ^)^ 2 ) + ^J^df- 2 + r 2 d<p 2 , (11) 
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where r + < r < oo and < (j) < 2n. Further defining 
the metric becomes, 



r 2 = rl(l + r -),t= r -^,x = r + ^ (12) 



2 2 

ds 2 = (1 + ^)rfx 2 + (1 + V + ^[-^ 2 + cosh 2 ^Q 2 ]. (13) 

We should note here that, in these coordinates, < r < oo. The Euclidean version 
of this space-time (t — > — iO — in) clearly resembles thermal AdS once we reinterpret 
the periodic coordinate x as Euclidean time. We record the metric here for later use: 

2 

ds 2 = (1 + f -)d^ + (1 + ^Y^dr 2 + r 2 [d0 2 + cos 2 M$. (14) 

Instabilities and decay rates: 

We now see from (13) and (7) that both these space-times have the same asymp- 
totic geometry. The boundary is time dependent and is given by dS 3 x S 1 . However, 
a notable difference is while the orbifold exists for any size of the asymptotic S 1 , the 
bubble appears only when this boundary circle has a maximum critical size. The size 
is given by the expression in (8). Boundary energy densities of the orbifold and the 
bubble are computed in [5]. They are given by 

P - hl = ~6^GV Pbubblc = -i^ (m+ I } - (15) 
Let us now consider the case where the size of the boundary circle is less than the 
critical value given in (8). In the bulk, we can have both the bubble or the orbifold 
geometry. However, in view of equation (15), we see that the orbifold will decay to 
the bubble of nothing by radiating away its energy [10, 11]. This is the analogue of 
Witten's decay of Kaluza-Klein vacuum to a bubble of nothing [6] 5 . 

As in the case of Kaluza-Klein decay, it is possible to find the bounce solution 
which mediates this decay. This was discussed in some detail in [11]. As analysed 
there, it is the Euclidean continuation of the smaller Schwarzschild black hole which 
acts as a bounce. The metric for the bounce is therefore 

ds 2 = (1 + £ - - 2 )d X 2 + (1 + £ - ^)-W 2 + r\d6 2 + costa 2 ). (16) 



5 When embedded in supcrsymmetric theory, one employs antipcriodic boundary condition of the 
fcrmions around the circle. This breaks supersymmetry completely. 
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Having identified the bounce, we can calculate the semiclassical decay rate from 
the orbifold to the bubble by evaluating the action difference between Euclideanised 
orbifold (14) and the bounce (16). Though individual actions diverge due to large 
volume, the difference remains finite once we require same asymptotic boundary con- 
ditions for x in (14) and for x in (16). This is obtained by setting 




4 + (17) 

Here (3^ is the period of x m (14) and f3 x is the period of x m (16); the expression of 
the later is given in (5). Now the difference in actions is given by 

AI -^bounce -^orbifold 

2x4 , A , fR „ , r _ rl% , r R 



16ttG/ 2 7o 

^3 ,HL/ 2 - r H 



J Q X dxJ r r 3 dr j ' dQ 2 3 - dx J q r 3 dr J dfl 2 3 } 



8G K 2rl + P 



(18) 



Here u 3 is the volume of unit three sphere. In getting the last expression we have 
made use of the boundary condition (17). 

Enhancement of decay rate due to string wrapping the circle: 

On generic grounds, we expect an enhancement of orbifold decay rate when a 
Nambu-Goto string wraps around the circle x- This is what we intend to compute in 
this subsection. We noted that the decay of the orbifold requires a bounce solution 
with a negative mode in its spectrum of small fluctuation. The small Euclidean 
AdS-Schwarzschild black hole that we analysed in the previous section has such a 
mode [24]. In what follows, we will assume that the presence of a string does not 
remove this non-conformal negative mode. With this assumption, it is now easy to 
see how the decay rate changes as we wrap a string with action 



S = T 



J d 2 tjd^f (19) 



along x direction of the orbifold and x direction of the bounce solution. Here, 7 is 
the induced metric on the string. We expect that the change in the decay rate will 
be proportional to the exponential of the action difference AS = bounce ~ ^orbifold- 
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This quantity can be easily computed as follows: 



AS 



bounce '-'orbifold 



T J d 2 eVdet7bubbie -T J c^^den 



orbifold 



= T 



/ dx 

JO Jr 



r+ 



R 



dr-T 



pit ( 

Jo Jo 



R 



dr 



2r 2 + + 1 2 ' 



(20) 



To get to the last step, we have used (17) and also made a large R approximation 6 . 
From the above expression we see that the decay rate increases with r + . This, in turn, 
implies that enhancement is larger for larger size bubble in the final state. We see from 
the above expression that AS increases with string tension T. However, for a string 
with large energy density, we would have to go beyond test string approximation. 

Decay via tachyon: 

In a recent paper [22] , Horowitz has argued that a black string can catalyse Wit- 
ten's decay process. When a black brane is wrapped around a compact circle, the 
circle size becomes a function of the position. For a suitable choice of brane, re- 
gardless of its asymptotic size, this circle can reach string scale at the brane horizon. 
In fact, tuning the charges of the branes, this circle size can be made to vary very 
slowly. For an anti-periodic boundary condition of the fermions around this circle, 
one expects a tachyonic mode to appear as the size shrinks to string scale. This mode 
may then induce a topology changing process by pinching off the circle at the horizon. 
This, in turn, creates a bubble. A concrete example of this is the D3 brane in ten 
dimension. For our purpose, instead of a flat space-time, let us consider N D3 branes 
filling up the boost orbifold R 1 ' 1 jZ. In the near horizon limit one gets AdS orbifold 
in five dimensions [11]. On the other hand, after the tachyon condensation, we have 
a bubble of the kind that we have been considering. Since this decay is catalysed 
by a string scale process, one would expect the rate to be much faster than the one 
through supergravity bounce. 

In the next section we discuss black holes, bubbles and orbifolds in GB theory. 
We study how the above features change as a function of the GB coupling. 

6 A similar computation was performed in [25] in the context of Witten's bubble 
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3 Gauss-Bonnet black holes 



We start by considering (n + 1) dimensional gravitational action in the presence of a 
negative cosmological constant A along with a GB term. 

R 



d n+1 x^j 



n+l 



- 2A + a{R 2 - 4R ab R ab + R abcd R abcd ) 



(21) 



This action possesses black hole solutions which we call GB black holes [20, 26-31] . 
In the above action, a is the GB coupling. As the higher derivative corrections are 
expected to appear from the a' corrections in underlying string theory, we will often 
refer to such corrections as a' corrections in this paper. The metric of these holes can 
be expressed as 



ds 1 



-V(r)dt 2 + ^- + r 2 dQ 2 n _ l , 



V(r) 



where V(r) is given by 



r 1 r 2 r 4a 4am 



V ^ = 1 + 2a-2a-l 1 -p + 



(22) 



(23) 



We first define various parameters that appear in the above equation. dfl 2 l _ 1 is the 
metric of a n — 1 dimensional sphere. I 2 is related to the cosmological constant as 
I 2 = —n(n — l)/(2/c n+ iA). Furthermore, we have defined a = (n — 2)(n — 3)aK n+ i, 
where K n+ i is the n + l dimensional gravitational constant. The parameter m in (23) 
is related to the energy of the configuration as 

(n - l)w ti _im 



M 



(24) 



where cu n _i is the volume of the n — 1 dimensional unit sphere. Asymptotically, the 
metric (23) goes to AdS space, since in this limit 



V{r) = 1 + 



1 1/46; 
L2q"~2q"v ~~P> 

We see from here that the metric is real if and only if 



a < I 2 1 4. 



(25) 



(26) 



In our discussion, we will always consider a satisfying the above bound. The metric 
(22) has a central singularity at r = 0. The zeros of V(r) correspond to the locations 
of the horizons. 



9 



In five dimension, for which n — 4, there is a single horizon at 

9 1 « 



-1 + 11 + 



4(m — a) 



/ 2 



(27) 



We note here that for a black hole to exist m > a. 

Thermodynamics of these black holes can be obtained via standard Euclidean 
action calculation. Such calculations were performed, for example, in [30] . Following 
these computations, the free energy and temperature can be written down as 



F = 



n-4 



K 



n+1 



n 



3)(r 2 + 2d) 



3)4(1-7?) 



6(n — l)dr\j 



+(n - 7)dr 2 + 2(n - 1)6? 



T 



[n 



AurAri + 2d) 



rl + 



n-4 



-a + 



n r] 



n-2 n-2 I 2 



The black hole entropy is given by 



/ 



T -i( dM 



47ru;„_ir 



n-l 



1 + 



12a 



n 



3r 



and the specific heat is 
dM 



C = 



dT 

47r(n — l)a; r , 



ir n- 3 ( r 2 + 2a) 2 [al 2 (n - 4) + r 2 (/ 2 (n - 2) + nr 2 )] 



(28) 



(29) 



K n+1 [ar 2 (6nr 2 - / 2 (n - 8)) + r 4 (nr 2 - (n - 2)l 2 ) - 2(n - 4)a 2 / 2 ] ' 

Many interesting features of the GB black holes, related to local and global stabilities, 
can be inferred from a detailed study of the thermodynamic quantities. In the rest of 
the section, we proceed to do so by considering the holes in five dimensions (n = 4). 
Let us first introduce two dimensionless quantities 



a = ^,and r = — . 



(31) 



We would like to express various thermodynamic quantities in terms of these dimen- 
sionless constants. The free energy given in (28) can be written as 

_ LU 3 l 2 

n 5 (f 2 + 2a) 

It then follows from (32), that within the range of allowed value of the coupling a 
(see (26)), F starts being positive at f = and changes sign only once as we increase 



f 6 + (18a - l)f 4 + 3af 2 - 6a< 



(32) 
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r. The number of extrema of the free energy, however, crucially depends on a. In 
particular, when a is in the region 

< a < ^, (33) 
F has three extrema. At these points, F takes non-zero positive values. However, for 

— < a < -, (34) 
36 ~ ~ 4' v ; 

F has no extremum for any non-zero f. It starts with a nonzero value at f = 0, then 
decreases monotonically and becomes negative at large f. Typical behaviour of the 
free energy as a function of f is shown Figure 1. We will refer back to this plot when 
we analyse the stability of these holes. 




Figure 1: Free energy as a function of x = r for different values of a. The thicker line is 
for at = 1/40 and the other one a = 1/32. 

For now, we turn our attention to the temperature of the black holes. It follows 
from (28) that the temperature is given by 7 

f + 2f 3 

T = 2vr/(f2 + 2a)- (36) 



7 In the limit I — ► oo (A = 0) this solution reduces to the asymptotically flat Gauss Bonnet black 
hole. The temperature is then given by 



2w(r 2 + + 2d) 



(35) 



For finite value of d, temperature begins with zero value at r + = and gradually increases for 
small r + . Finally it reaches a maximum value at r + = \[2a. and then again goes towards zero at 
large r + . Since the temperature has a maximum, above this critical value there is no black hole 
solution. At any temperature below there are two black hole solutions, small and large. The small 
black hole has positive specific heat and is locally stable. The larger one is unstable due to its 
negative specific heat. This is to be contrasted with the Schwarzschild black hole solution without 
R 2 correction, where we have only one unstable solution existing at all temperatures. 



11 



At f = 0, temperature starts out from zero and, regardless of the value of a, it 
increases for small f. However, at larger f, the number of extrema depends on a. In 
the region given in (33), there are two of these extrema. Both of these disappear as 
we increase a to region (34). A plot of the temperature as a function of f is shown 
in Figure 2. 

T 

3.5 

3 

2.5 

2 

1.5 

1 

0.5 

0.25 0.5 0.75 1 1.25 1.5 1.75 " 

Figure 2: Temperature as a function of x = f for different values of a. The thicker line is 
for a = 1/44 and the other one a = 1/30. 

To examine the phase structure of these black holes, it is instructive to consider 
the behaviour of the free energy as a function of temperature (for different values of 
a). From (32) and (36), it is possible to construct the temperature dependence of the 
free energy. However, the analytical expression is not very illuminating. Therefore, 
we plot the nature of the free energy as a function of temperature in Figure 3. This 
plot is for two different values of a belonging to the two different regions given in 
(33) and (34). Note that, as we increase a from region (33) to region (34), nature of 
F changes at a critical value a = a c = 1/36. We therefore study these two regions 
separately. 

Phase structure for a < a c : 

When a < a c , the free energy is shown by the thicker line in Figure 3. At low 
temperature, it has only one branch (shown as branch I in the figure). However, 
when the temperature is increased beyond a certain value (which we call Ti), two 
new branches appear (II and III). One of these two branches (II) meets branch I at 
a temperature beyond, say T 3 , and they both disappear. On the other hand, branch 
III continues to decrease rapidly, cuts branch I at temperature, say T2, and becomes 
negative at a temperature which we will call T c in the future. While computing 
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Figure 3: Free energy as a function of temperature. The thicker one is for a = 1/50 while 
the other one is for a = 1/30. 

specific heat using (30), we find that it is positive for branch I, and III. Therefore, 
these phases correspond to stable black holes. They, however, differ in their sizes; 
branch I represents smaller sized black holes than that of branch III. Going back now 
to branch II, we find that the specific heat is negative. We, therefore, conclude that 
branch II represents an unstable phase of the black hole. 

The above picture is similar to that of the van der Waals gas. In particular, the 
Gibbs free energy of van der Waals gas, for an isotherm, behaves in a similar manner 
as we vary pressure. A thermodynamic equilibrium state is reached by minimising 
the Gibbs free energy. Likewise, in our case, equilibrium state would correspond to 
branch I of the free energy all the way up to temperature T 2 and then branch III from 
temperature T 2 and above. The free energy curve then remains concave as expected 
for a thermodynamical system. We, however, note that since there is a discontinuity 
of dF/dT at T = T 2 , one has a first order phase transition at T 2 . Two black hole 
phases would differ from each other at this point by a discontinuous change in their 
entropies. We will call these as the first Hawking-Page (HP1) transition for reasons 
that will be obvious later. 

This phase structure can be nicely described by constructing a Landau function 
around the critical point. By identifying the dimensionless quantity f as an order 
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parameter, we can construct a function <3>(T, r) as 8 

U) l 2 

$(T, f) = — (3f 4 - 47rlTf 3 + 3f 2 - 2AnalTf + 3a). (37) 

At the saddle point of this function, that is when ^ = 0, we get back the expres- 
sion of the temperature given in (36). If we then substitute back the expression of 
temperature in to (37), $(f) reduces to the free energy given in (32). As can be 
seen from Figure 4, for temperature T < T 2 , $(T, f) has only one global minimum. 
This corresponds to the small black hole phase. However, at T = T 2 , appearance 
of two degenerate minima suggests a coexistence of small and big black hole phases. 
Finally for temperature beyond T 2 , only the big black holes phase remains (as this 
phase minimizes the Landau function). Clearly, there is a discrete change in the order 
parameter f at T — T 2 . This is what we expect for a first order phase transition. 

Phi 




-0 . 02 [ 

Figure 4: Landau function $ as a function of order parameter x = f for different temper- 
atures. We have taken a = 1/50. 



Phase structure for a > a c : 

For a > a c , the free energy curve is shown by the thin line in Figure 3. Unlike the 
previous case, free energy and its derivatives do not show any discontinuity. Therefore, 
there is no HP1 transition. Only a single black hole phase is found to exist at any 
temperature. 

Global phase structure of GB black holes: 

As discussed earlier, for black holes with a — 0, there is a crossover from AdS to 
AdS black holes at a critical temperature ^7. What happens to this transition as we 

8 To construct the Landau function, we employ a method similar to the one discussed it [32,33]. 
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turn on a! In this situation, we note that we still have two geometries to consider. 
First one is again a thermal AdS with metric being the Euclidean continuation of 
(22). The function V(r) is given in (25). We identify this thermal AdS space, having 
a dependent effective cosmo logical constant, with zero free energy. Now, from Figure 
3. we see that above a critical temperature, the free energy of the GB black hole 
becomes negative, making it more stable compared to the effective AdS geometry. 
We identify this as a HP2 point. This crossover temperature can be computed as a 
power series in a and is given by 

We notice here that the GB correction reduces the transition temperature. Similar 
phenomenon was noticed earlier in many AdS-gravity theories with higher curvature 
terms [34-37]. The global phase structure is shown in Figure 5. 

F 

0.4 
0.2 



-0.2 
-0 . 4 
-0 . 6 
-0 . 8 

-1 

Figure 5: Global phase structure of GB black holes. For temperature T < T c , AdS lowers 
the free energy, on the other hand for T > T c the black hole phase is preferred. This plot 
is for a = 1/30. 

To this end, we would like to point out that the above picture of GB black holes 
is quite similar to that of the five dimensional charged AdS black holes [38,39]. The 
stability properties of the charged black holes depend on whether we are considering 
fixed potential ensemble or fixed charge ensemble. For the case of fixed charge en- 
semble, various phases of black holes resemble that of the thick line in Figure 3. As in 
our case, small black holes and large black holes are separated by a first order phase 
transition point. However, a major difference is that for the charged black holes, in 
fixed charge ensemble, thermal AdS is not a solution. Consequently, these holes are 
globally stable. This is unlike GB black holes, where there is a HP2 transition. Below 
HP2 temperature, they are unstable. 



T_C 



0.5 
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4 Bubbles in GB theory and instabilities 



In this section, we turn our attention to bubble spacetime in GB theory. After 
constructing these bubbles, we compute their energy densities. We find that due to 
the presence of these bubbles with same asymptotic structure as of the AdS orbifolds 
in GB theory, the orbifolds become unstable. The nature of these decays depend 
on the value of a. In what follows, we will be mainly interested in studying these 
instabilities as a function of a. We will also highlight the differences that occur when 
we compare the present situation with the one with a = 0. 

The bubbles in the GB theory can be constructed by analytically continuing the 
coordinates of the GB black holes as t — > i X , 6 tt/2 + it. Here theta parametrises 
one of the angles of d£l n _ x in (22). The solutions then takes the form 

dv^ 

ds 2 = V(r)d X 2 + - r 2 dr 2 + r 2 cosh 2 l rdQ 2 n _ 3 , (39) 

where V(r) is given in (23). Above metric is nonsingular in the region r > r + if X 
has a periodicity 

Anfl(f 2 + 2a) 

*= ( n _2)(f2 + 2 E |a + ^)- ( } 

In the following we will use dimensionless quantity Ax = A X /l to parametrise the 
circle. We first note that at asymptotically large distance, the metric reduces to 
dS n -i x S* 1 where S 1 corresponds to the X circle. More precisely, up to a conformal 
scaling by L 2 /r 2 , the boundary metric becomes 

ds 2 = d X 2 + L 2 (-dr 2 + cosh 2 rdn 2 n _ 2 ). (41) 

In the above equation, we have defined 

L = V2& 



/ 2 > 

At r = r + , the proper radius of this circle collapses at V(r) = 0. However, the 
n — 2 sphere approaches a finite size r\ coshV. Therefore (39) represents a bubble of 
nothing in GB theory with size r+. In the rest of this section, we will mostly focus 
ourselves on the bubbles in five dimensions. 

Many of the features of these bubbles can be understood from the behaviour of 
the periodicity Ax as a function of a and f. Firstly, as can be seen from Figure 6., 
for any non-zero a, there exists a bubble for any size of the X circle. This is very 
much unlike the case when a = where there is a critical radius above which the 
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bubbles are no longer present. Secondly, for a given a in the range (33), and for a 
given periodicity of Xi there can be at most three bubbles of varied sizes. This can 
be seen from the solid line in Figure 6. However, as we increase a above 1/36 and go 
to the range (34), for fixed Ax, we get a single bubble of fixed size. This is shown by 
the dashed line in Figure 6. 

chi 

5\ 



Figure 6: Plot of Ax as a function x = r. The dashed line is for a = 1/34. The solid line 
is for a = 1/50 and the dotted line is for a = 0. 



Continuing our discussions in five dimensions, let us note that as in the case of 
AdS gravity, in GB theory, we also have another geometry with the same asymptotic 
metric. These are the AdS orbifolds (13) with the AdS curvature / replaced by L as 
in (42). 

The asymptotic boundary of the GB bubble spacetime is dS 3 x S 1 . One would 
then expect, by AdS/CFT correspondence, that some deformation of M = 4 SU(N) 
Yang-Mills theory at large but finite A should be dual to the bubble. Clearly, similar 
to the AdS bubbles, for the GB bubbles, the CFT lives on a time dependent space. 
The boundary stress tensor can be computed from the bulk stress tensor. These bulk 
stress tensor components, for GB bubbles, can easily be obtained and are given by 

1 



'TX 

X 



T 



rpq 



m, 



K 5 L 3 

3 



m, 



K 5 L 3 

1 

K 5 L 3 



m, 



m, 



(43) 
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where we have parametrised dVl 2 by the coordinates 9 and 0. We note here that the 
components of the stress-tensor computed with respect to the orbifold in GB theory 9 . 
Positive sign of implies that the solution has negative mass. 

As of now, we have learned that, in GB theory, bubble space exists for any value 
of m. This, in turn, means that we can have a bubble of any size. Furthermore, we 
note that the bubble and AdS orbifold have asymptotically the same metric, namely 
dS s x S* 1 . The energy density of the bubble spacetime is however less than that of the 
orbifold. We may, therefore, conclude that AdS orbifold is an unstable background in 
GB theory. It will always decay to bubble by radiating away its energy. As a result, 
bubble of any size will be produced. Once it is produced, due to the time-dependent 
nature of the metric, the radius of the bubble will increase exponentially with time. 

It is easy to compute the decay rate from the orbifold to bubble. This can be 
done, as before, by identifying the bounce solution. 

5 Matrix Model : Some numerical computation 

In the previous sections we have analysed the phase structure of the gravitational 
theory in the presence of a higher derivative correction. In the light of the AdS/CFT 
correspondence we would now like to analyse these phases from the gauge theory on 
the boundary. Specifically, we will study the thermal aspects of the SU(N) gauge 
theory on the boundary of AdS 5 which is S* 3 x S 1 . The effective theory on the 
boundary can be described by an unitary matrix model: 



where U = Pexp(i Jjf A(r)dr) and A(r) is the zero mode of A on S* 3 . This is the 
lightest mode, and the effective action is obtained by integrating out all the massive 
modes. In general S e ff(U) is a polynomial in the traces of U and its powers that are 
allowed by the Zjy symmetry. The coefficients of these terms depend on the 'tHooft 
coupling A (that is related to a by, a V2X = I 2 from the AdS/CFT correspondence) 
and temperature T. 

These coefficients have been worked out in the weak coupling expansion in [41]. 
In our analysis we will restrict ourselves to the first two terms, 



9 One easy way to compute these expressions is to first calculate various components of stress- 
energy tensor for GB black holes. The formalism to compute for the GB black holes is given 
in [40]. One can then make a proper analytical continuation to get stress tensor for the GB bubbles. 




(44) 



S eff (U) = atrC/trC/ 1 + -^-(trUtrU^ 2 , 



(45) 
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where, a and b are functions of temperature T and A. An order parameter charac- 
terising the deconfined phase of the gauge theory is given by the expectation value of 
the Polyakov loop l/N(trU). 

The saddle point equations for this effective theory are given by, 

1 

ap + 2bp s = p < p < - 

\<P<^ (46) 



4(1 -P) 2 

where p 2 = (l/N 2 )trUtrW . As mentioned this matrix model contains only the first 
two terms of the effective gauge theory in the weak coupling expansion. It was shown 
in [19] that with this truncation, in the large N limit one can reproduce the same 
thermodynamic features as of the bulk black hole thermodynamics near the critical 
points. This model thus appears to fall in the same universality class as that of the 
boundary effective gauge theory in the strong coupling limit. 

The phase structure in the bulk theory that we have discussed in Section 3 con- 
tains various distinct qualitative features depending on the value of the correction 
parameter a . For nonzero a , the phase diagram is modified in the regime where 
r + is small compared to \fd . However as long as r + (the solutions corresponding 
to the black holes at a particular temperature) are greater than a', the phase dia- 
gram is qualitatively the same as that of the bulk theory without higher derivative 
corrections. There are two possibilities. 

• We can ignore this small black hole solution in the supergravity approximation, 
so that we only concentrate on solutions r + > \fa' ' . In this domain, it makes 
sense to compare the bulk physics with that of the boundary (a, b) matrix model 
discussed in the earlier paragraph. 

• If we include the small black hole solution in the supergravity approximation, 
then in order to reproduce the bulk phases, the boundary matrix model needs 
to be modified. In the next section we will propose a matrix model potential 
that captures the bulk physics including the solution r + which is less than \fa! \ 

The following part of this section is devoted to the study of the (a, b) model 
numerically, incorporating the corrections due to the finite 'tHooft coupling A. In 
this analysis we will take N — > oo. We first work in the limit A — > oo and then by 
taking A large but finite. The main aim is to compute a and b as functions of T and 
of A (to the first order in 1/ y/X) . This will be done by comparing the matrix model 
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potential with the action on the gravity side with a' corrections. In this paper we 
are considering the corrections due to the Gauss-Bonnet term. As mentioned before 
R 2 corrections are not known to occur in the supergravity limit of type IIB theory. 
However the following analysis is a fruitful exercise that can easily be adapted for the 
R 4 terms that arise in this theory. 

The comparison between matrix model potential and the action of the bulk theory 
is valid as long as we can neglect the string loop corrections. The corresponding 
temperature at which the supergravity description breaks down is identified as the 
Gross- Witten transition point in [19,43] on the matrix model side. 

Let To be the temperature at which the black hole nucleation starts. For T > To, 
it is well known that for the gravity theory without a' corrections, one gets two 
solutions for the black hole. The small black hole is unstable and the larger one 
stable. The larger one undergoes a Hawking-Page transition at T = T c . It was shown 
by Witten [3] that thermal AdS§ corresponds to the confined phase of the large N 
gauge theory on the boundary. A natural order parameter that characterises the 
deconfmed phase is the Polyakov loop. In this matrix model it is p. 

Let us study the case without a' corrections first. This corresponds to the A — > oo 
limit. Consider T > T , for which we have 10 



2ap\ 2 + 2bp{ 2 + log(l - p 1)2 ) + / = -J 1>2 , (47) 

where the are the actions for the large and small black-holes respectively and p 12 
are the corresponding solutions in the matrix model. The constant / = log(2) — 1/2 

10 We will set I, W3, k 5 to 1 in the numerical computations. 
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Figure 8: The main figure is the plot of left (dashed lines) and the right hand side 
(solid line) of the saddle point equations. T is the point where the two roots of (46) 
merge. T c is the curve corresponding to the Hawking-Page transition temperature. 
The inserts A and B are the potentials corresponding to these temperatures 



is added to make the potential from (46) continuous at p = 1/2. Since the values 
of pi j2 are those at the extremum of the left hand side of (47), we have two more 
equations that are given by (46). 

For a given temperature, T, are known from the gravity side, so the problem 
now is to solve the above equations for a(T), b(T) and p lj2 . We do this numerically. 
The solutions are plotted in Figure (7). Note that a(T) and b(T) increases monoton- 
ically. As a consequence of this the dashed line in Figure (8) that represents the left 
hand side of (46) moves towards T c as the temperature is increased thus generating 
two solutions for p corresponding to the small and the big black holes in between. 
p = corresponding to the thermal AdS$ is always a solution. 

Having known the variations of a(T, 0) and b(T, 0) with respect to the temperature 
we now incorporate the a' corrections to to get the first order dependence on 
1 / We have 

a(T,l/V\) = a(T,0) + -)=^L-| 1/v ^ =0 +O(l/A 3 / 2 ), (48) 



21 



Figure 9: Plots of (A) da(T)/d(l/y/\) and (B) db(T)/d(l/y/\) 



The first order variations of equation (47) with respect to 1/ y/X gives, 



2-^P? 2 + 2^^PU = -VXSh 2 (T). (49) 
In the above equations, 5I\^{T) are given by, 



5I lt2 {T) = dp(6F lt2 ) 

j :(3rf 2 + 24r 1 2 2 + 9). (50) 



'2A 

Where F is given by eqn(32). From these we get the values of q^J^ an d d(\j^}x) 
shown in Figure (9). 

At this point some comments about the sign of b are in order. In the limit when 
the t'Hooft coupling A goes to infinity numerical computations show that b indeed is 
positive. However as we move from A — > oo to finite A (that is obtained by including a' 
corrections in the bulk) we see that at any particular temperature db(T)/d{l/ \J\\)) 
is negative. Though this numerical calculation shows that b decreases from a positive 
value as we move down towards weak coupling, it is not clear whether the sign of 
b will turn out to be negative or positive at weak coupling. In case it turns out to 
be positive, we presume that the results for the gravitational side should correspond 
qualitatively to those of the weakly coupled gauge theory. 
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The above analysis shows that the behavior of the coefficients as functions of 
temperature and A are indeed the ones that we expect from the phases of the bulk 
theory as long as we concentrate on the black hole solutions with r + > \fd . The 
expansions are carried out about A — > oo as it was argued in [19] that the effective 
theory that is computed in the weak coupling falls in the same universality class as the 
one in the strong coupling limit. The addition of higher derivative term in the bulk 
does give information about l/y/\ corrections, however this (a, b) model is unable to 
capture the phases including the small black hole solution. In the following section 
we will analyse this issue, in detail, by proposing another model which qualitatively 
reproduces various bulk phases of section 3. 



6 A modified Matrix model 

In this section we propose a modified (toy) matrix model which incorporates some of 
the additional qualitative features on the gravity side arised due to the GB term. We 
find, the minimal action that would reproduce these features needs to be quartic in 
p 2 and can be given by 

i _oi 

S(p 2 ) = 2[A 4P 8 - A 3P e + A 2P ' + (^-^)P 2 ] , (51) 

where A^s are the parameters, which depend on the temperature as well as on the 
coupling constant. In the limit where the A4 and A 3 vanish we get the (a, b) model [19]. 

The equations of motion ensuing from the action in (51) are given as follows. We 
write 



F (P) = = t 8 ^P 6 - 6A *P A + 4 ^P 2 + (1 - 2A X )\. (52) 

Then the equations in two different regions are 

pF(p) =p , 0<p<±, 

= WT) - §<P<1- (53) 
The potentials that follows from the above action is given by 

V(p) = -A 4 p 8 + A 3 p 6 - A 2 p 4 + A lP 2 , 0<p<i , 

= -A 4 p 8 + A 3 p 6 -A 2 p i + (A 1 - 1 -)p 2 -hog[2(l-p)] + ± , \<P<1 

(54) 
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Let us analyze the solutions of equations of motion given by (53). The fact that 
there are four parameters instead of two has made the analysis technically more 
involved than (a, b) model [19]. For various ranges of parameters the model shows 
different qualitative behaviour. As we will see we need to impose necessary restrictions 
on the parameters so that the model reproduces the features that we found on the 
gravity side. Before analyzing the solutions one comment is in order. In the following 
we will find the analog of small stable black hole appearing as a minimum of the 
potential but it always comes with an additional maximum of the potential. We do 
not have on the bulk side a solution corresponding to this maximum. We interpret 
this solution as a possible decay mode of the small stable black hole which may be 
due to some stringy mechanism. 

The behaviours of the solutions are encoded in the polynomial F(p) given in (52). 
We begin with the coefficient of the lowest order term A\. From (51) we see in order 
to make p = tachyon-free we need < F(0) < 1 i.e. < A± < 1/2. Once that is 
imposed we consider the next coefficient A 2 . As we see on the bulk side our action 
should admit ( in one phase) 3 solutions that correspond to a small black hole, an 
intermediate black hole and a big black hole. A necessary condition for the existence 
of three solutions is A 2 > 0. Though we get this constraint from a different argument 
it agrees with [19]. Thus our model at a vanishing limit of higher coefficients reduces 
to (a, b) model. 

Restrictions on the higher coefficients are slightly more cumbersome and depends 
on the positions of the turning points of the polynomial F. In that context it is useful 
to consider the quadratic polynomial in p 2 : f(p 2 ) = (1/ p)-§-F(p). This is given by 
f(x) = A8A 4 x 2 — 2AA 3 x + 8A 2 . The zeroes of / determine the non-trivial turning 
points of F. In terms of this polynomial / the two different ranges of a correspond 
to the following constraints: 

• a > 1/36: For /(l) = A8A4 - 24A 3 + 8A 2 < there is one turning point at 
some < p_ < 1. With parameters in this range we can have either two solutions 
(one maximum and one minimum of potential) or no solution. There is no way 
we can obtain three solutions in this phase. Moreover, from the restriction on the 
parameters it is clear that the range of parameters is not continuously connected 
with the corresponding range where (a, b) -mo del is valid (i.e. A^ = A 3 = 0). We 
identify this phase with the range of a which corresponds to a > 1/36 on the gravity 
side. However, that is not sufficient to ensure that there is always one minimum that 
correspond to the single black hole on the bulk side. For that we need to impose a 
further restriction on the coefficients such that, the turning point satisfies p_ < 1/2 
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and F(p-) > 0. Then we always get a maximum for p < 1/2 ( that is in the region 
with no cut) and a minimum of the potential. As the parameter varies the position 
of this minimum changes from the p < 1/2 region to the p > 1/2 region. So this 
phase corresponds to the restrictions: /(l) = 48 A 4 — 2AA 3 + 8A 2 < 0, p_ < 1/2 and 
F(p_) > 0. 

• a < 1/36: Again we look for turning points in the range < p < 1. For 
/(l) = 48 A A — 2AA 3 + 8A 2 > either we get two turning points which we call p_ 
and p + (p_ < p + ) or none of them. That gives rise to 3 possibilities: the number of 
solutions could be 4 (consists of two maxima and two minima), or 2 (consists of one 
maximum and one minimum)or 0. This phase is continuously connected with that of 
the (a,b) model and we identify this phase with the range of a given by a < 1/36 
on the gravity side. The more detailed structure of the solutions depends on the 
position of the turning points p_ and p + and the values of the polynomial F(p) at 
the turning points. Let us first consider F(p_) > with < p_ < 1/2. There are 
two possibilities: (i) If _F(l/2) < we have two solutions, one maximum and other 
minimum in p < 1/2 range. The minimum corresponds to the stable small black hole. 
We may or may not have two more solutions in the range p > 1/2. If we have two 
solutions they would correspond to intermediate and big black hole, (ii) If F(l/2) > 
we have one solution (maximum) in p < 1/2 and the other (minimum) in p > 1/2. 
This minimum corresponds to the big black hole. The remaining possibilities are (iii) 
p_ < 1/2, F(p_) < and (iv) p_ > 1/2, F(l/2) < 0. In both of these cases there 
is no solution in the range < p < 1/2. Finally if there is no turning point and 
F(l/2) < there will be no solution in the range < p < 1/2. Since the analysis on 
the bulk side then requires that there is a solution for 1/2 < p < 1 we need to impose 
the following constraint, namely, there should exist some value of p, p„ < p < 1 
such that 4p (l — p )F(p ) > 1. That will give one maximum and one minimum in 
the range 1/2 < p < 1 that corresponds to the small and the big black hole. 

Thus we see for both the phases we need additional restrictions which shows there 
are regions of parameters that does not agree with the features of gravity phase. This 
suggests the fact that in the strongly coupled gauge theory there are restrictions on 
various parameters. It may be interesting to calculate these parameters from field 
theory set up (for weakly coupled gauge theory) and compare the values with the 
restrictions obtained above. 

In order to discuss the variation of potential with parameters it is useful to give 
a graphical description. We have four parameters, so for the sake of graphical de- 
scription we restrict number of parameters to 2. We consider only the phases that 
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corresponds to the a < 1/36. We take fixed values of A 1 and A 4 and study the 
features with the variation of two other parameters. We have chosen the values to be 
Ai = .025 and A 4 = 2.083. We have given plot of the potential against p in Fig. 10 and 
11. In order to make the extrema explicit we choose different scales for the potential 
for two different ranges of p, namely, < p < 1/2 and 1/2 < p < 1. The values of 
A 2 and A 3 are decreasing from the curve in bottom to that in top in Fig. 10 and from 
the curve on top to that in bottom in Fig. 11. There is always one minimum at p = 
where the potential vanishes. As we will see in Fig. 12 we need to choose values of A 2 
and A% restricted within a particular region outside which the features that we get 
from the bulk will be absent. In the following, we give V(p) vs. p plots for different 
values of A 2 and A 3 : 

• (A 2 ,A 3 ) = (.45,2): Here we get two solutions: one maximum and one local 
minimum in the range p < 1/2 (i.e where there is no cut) and no solution at 
p > 1/2. We identify the minimum with the small stable black hole. This 
corresponds to low temperature behaviour of GB black hole where we get only 
one small black hole solution. 

• (^2,^3) = (.4,2): Here we get four solutions: In addition to the above maxi- 
mum and minimum in the range p < 1/2 we get one more local maximum and 
one more local minimum appearing in the range 1/2 < p < 1. These latter 
maximum and minimum can be identified with the intermediate black hole and 
the big black hole. We identify this with the nucleation of the big black hole 
and intermediate unstable black hole in the gravity picture. 

• (A 2 , A 3 ) = (.385, 1.9375): For further decrease of the parameters, the heights of 
the local minimum in the range > p > 1/2 (Fig. 10)increases and the height 
of the local minimum in the range 1/2 > p > 1 decreases(Fig. 11) . At this 
value of the parameters the heights of the two minima become equal. We can 
identify this point with a transition from small black hole to big black hole on 
the gravity side which is termed as HP1 transition. 

• (A 2 ,A 3 ) = (.38485, 1.93688): (Due to close proximity this plot appears on the 
top of the earlier plot and not distinguishable in the present scale.) The height 
of the minimum in the range p > 1/2 becomes zero and thus equal to the 
potential at p = 0. On the gravity side this corresponds to energy of big black 
hole reaching zero and becoming equal to that of thermal AdS triggering HP2 
transition. 
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• (A 2 ,A 3 ) = (.25, 1.5): Here we get two solutions because in the region p < 1/2 
the local minimum and local maximum is on the verge of disappearing. However, 
the two solutions in the range p > 1/2 will remain with the height of the 
minimum in p > 1/2 keeps on decreasing. This corresponds to the point beyond 
which the small black hole on the gravity side disappears. 

• (A 2 ,A 3 ) = (.248, 1.25): As we decrease A 2 and A 3 further, the solutions in the 
range p < 1/2 cease to exist (Fig. 10). The minimum in the range p > 1/2 (Fig. 
11) becomes more and more deeper. This is in keeping with the fact that, at 
high temperature on the gravity side the only stable configuration remains is 
the big black hole. 

As we see from the above analysis the coefficients decrease with temperature, 
unlike the behaviour of the coefficients in the (a, 6)-model. This can be interpreted 
as the temperature gradient of the coefficients at the first order of inverse 't Hooft 
coupling has a negative sign relative to that at the zeroeth order. At this range, 
where appreciable 1/A correction is taken into account, the contribution at first order 
dominates over that at zeroeth order. 

(A2, A3) 




0.1 0.2 0.3 0.4 0.5 

Figure 10: Potential as function of p for the range < p < 1/2 for increasing values of 
A 2 and ^3. The different values of (A 2 , A3) are given above. The plots associated with 
(.38485, 1.93688)and (.385, 1.9375) are not distinct in this scale. 
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( A2 , A3 ) 
(.248 , 1.25 ) 




Figure 11: Potential as function of p for the range 1/2 < p < 1 for increasing values of 
A 2 and A3. The values of (A2, A3) used in the plots are given above. The plots associated 
with (.38485 , 1.93688) and (.385 , 1.9375) are not distinct in this scale. 

Here we give a graphical presentation of the behaviour of the solutions using a 
parametric plot of different critical points in the A 2 -A 3 plane in Fig. 12 keeping Ai 
and A^ fixed as above. As we vary the parameters we encounter a curve IV in the 
A 2 -A 3 plane, above which the saddle point associated with the small black hole has 
energy negative. From the analysis of black holes on the gravity side, it follows that 
the small black hole energy is always greater than thermal AdS ensuring the stability 
of the latter. So, in what follows, we restrict ourselves to the region below curve IV. 

In the region bounded by IV, III and I, there are three saddle points. One is p = 
which corresponds to the thermal AdS. There are two more saddle points: a local 
maximum at p — p\ and a local minimum at p = p 2 . The latter corresponds to the 
small black hole that we obtain on the gravity side. There is no solution analogous 
to p\ in the gravity side. In the region bounded by II, III, IV and I, there appears 
two more saddle points. One of them p = p 3 is a local maximum and the other one 
p = Pa is a local minimum. They correspond to the intermediate, and the stable 
big black hole respectively. In the region on the left hand side of curve I, the saddle 
points p = pi, P2 cease to exist. In the region above the curve IV, as we have already 
mentioned, the potential of the saddle point p = p\ becomes negative showing the 
energy of the associated small black hole on the gravity side becomes less than that 



28 




1.25 1.5 1.75 2 2.25 2.5 2.75 3 

Figure 12: Parametric plots of different critical points in the A2-A3 plane. We choose 
Ai = 0.25 and = 2.083. In the region which is below or on the left side of curve I the 
saddle points p± cease to exist. In the region above curve IV the potential of p + vanishes. 
Curve II corresponds to HP transition and on curve III the saddle points p\ , P2 merge. 

of thermal AdS. 

Similarly, the thermal history (for this choice of parameter) can be obtained from 
Fig. 12 as follows. As we mentioned earlier, A 2 and A 3 will decrease with temperature 
along the curve C. As we follow the curve C from right to left, we find the p = 0, pi, p 2 
are the solutions on the right of curve III. As we cross curve III, we encounter two 
additional saddle points p = P3, p^. Crossing the curve II corresponds to the Hawking- 
Page transition. As we cross curve I, the saddle point corresponds to the small black 
hole disappears. 

Like the general case, here also as we see in the region bounded by the curve 
I, along with a local minimum ( at p — pi) we always obtain a local maximum ( at 
p = P2). We interpret this maximum, as we said earlier, as a bounce solution through 
which the small stable black hole decays. It will be interesting to understand this 
instability on the gravity side. 
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7 Discussion 



In this paper, we have discussed phase transition of asymptotically AdS black hole 
solutions in presence of Gauss-Bonnet term. As long as a, strength of the coupling to 
GB term remains above certain critical value a c , one gets a single black hole phase at 
any temperature. However, as the coupling comes down below the critical value two 
additional black holes appear. We called them small and intermediate black holes. 
The intermediate black hole is found to have negative specific heat. It turns out that 
this small stable black hole is a local minimum below a critical temperature. Beyond 
this temperature small black hole disappears. We have studied the associated phase 
diagram and find that the phase structure resembles that of van der Wall's gas. In 
addition to the the standard Hawking-Page transition, we have identified one more 
phase transition where the two branches of the phase diagram meet. We find the 
specific heat diverges at this new critical point. 

From a different perspective, the Euclidean version of the black hole solution has 
been interpreted as the bounce mediating the decay of AdS orbifold into the bubbles 
of nothing. It was found earlier that for the radius of circle above a critical value the 
AdS orbifold is stable while below that bound the AdS orbifold decays. After adding 
the higher derivative terms in the action, we find that there is a bubble solution that 
exists for all values of the radius and has an energy density lower than that of AdS 
orbifold. This suggests that the AdS orbifold is unstable for any radius. One of the 
stringy feature that we have addressed in this paper is the modification of decay rate 
in presence of a background string wrapped on the circle and we find the decay rate 
gets enhanced. It would be interesting to find tachyonic decay ensuing from wrapped 
black brane in this set up. 

Five dimensional theory of gravity usually corresponds to some gauge theory on 
the boundary and the analysis on the gravity side has natural implications about 
the gauge theory. In absence of Gauss-Bonnet term, the gravity theory on Euclidean 
AdS ( along with S 5 ) is known to be dual to be pure N=4 SYM on a three sphere at 
finite temperature and the phase diagram associated with the gravity theory captures 
thermal history of N=4 SYM on S 3 . In a similar spirit, we expect, dual of this five 
dimensional gravity theory in presence of Gauss-Bonnet term is some deformation of 
the above gauge theory and the phase diagram captures its thermal history. In [19], 
qualitative features of N = 4 SYM on S 3 was studied from the perspective of a 
matrix model. This model is phenomenological in nature and is charecterised by two 
parameters (a, b). On generic ground, one expects these parameters to be A and T 
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dependent. Appealing to the universal nature of this model near the critical points, 
we find out A dependence of (a, b). This is done by mapping the bulk a' correction 
to the boundary. This method can easily be used to find similar 1/A dependence 
of matrix model coefficients in the case of other higher derivative corrections of the 
gravity action, such as, R 4 term in IIB theory. 

We have also proposed a modified matrix model that captures the qualitative 
features of the phase diagram of the bulk theory. Unlike (a, b) matrix model this 
model is non-universal and the phase diagram is reproduced only in a selected region 
of the parameter space. In addition the temperature dependence of the coefficients 
turn out to be different from usual linear increasing function. We also find there is a 
bounce through which the small black hole can decay. It will be interesting to identify 
this instability on the gravity side. We hope to return with some of these issues in 
future. 
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